Voltage-controlled spin injection with an endohedral fullerene Co@Ceo dimer 
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Spin-dependent transport through an endohedral fullerene Co@C60 dimer with gold electrodes is 
explored theoretically using density functional and extended Hiickel theory. Density of states spin 
polarizations up to 95%, due to spin-splitting of Co 3d orbitals, are found by varying the gate and/or 
bias voltage. The current-voltage characteristics and strong (up to 100%) spin polarization of the 
current indicate that the device can be utilized for highly efficient spin injection into nonmagnetic 
conductors. This finding opens the way to the realization of electrostatically tuned spintronic nano 
devices less than 2 nanometers in size, without ferromagnetic electrodes. 



Carbon-based nanostructures such as fullerenes, car- 
bon nanotubes (CNTs), and graphene, are promising 
candidates for spintronic applications because of their 
weak spin-orbit coupling and hypcrfine interaction which 
lead to long spin coherence lengths [1, 2]. In particu- 
lar, the fullerene Cgo molecule is an interesting carbon 
nanostructure that can be used as a molecular bridge in 
magnetic tunnel junctions due to its remarkable struc- 
tural stability and electronic properties [3, 4]. 

One way to generate a spin-polarized charge current 
is to encapsulate magnetic atoms or magnetic nanowircs 
in fullerenes or CNTs [5-20]. For example, ab initio cal- 
culations showed that Co atoms in CNTs can provide 
strong spin polarization and considerable magnetic mo- 
ments [5]. The electrical and magnetic properties of Ceo- 
Co nanocomposites have also been studied [7, 11, 12] 
and tunnel magnetoresistance ratios up to about 30% 
at low bias voltages were reported, indicating significant 
spin polarizations [7]. Recent density functional theory 
(DFT) studies of encapsulated Co atoms in Ceo molecules 
[18] found hybridization between Co and Ceo orbitals and 
the most stable structure of Co@Ceo to have the Co atom 
on top of a hexagonal face. Lu et al. [19] showed that, 
in the most stable structure of Gd@C6o, the Gd ion is 
over a hexagonal ring of the Cgo molecule and the Gd 
atomic orbitals to hybridize with the Cgo molecular or- 
bitals. Moreover, strong hybridization between the Gd 
5d and 6s and carbon orbitals in Gd@Cgo has been ob- 
served, both theoretically and experimentally [20]. 

Magnetic atoms encapsulated in carbon nanocages, 
like the Cgo molecule, are effectively protected against 
environmental effects such as oxidization, stabilizing the 
encapsulated magnetic atoms for potential applications 
for spin injection in nanoscale devices. Accordingly, one 
can propose magnetic nano junctions to produce high 
spin-polarized currents with long spin coherence lengths 
[21]. In this paper we explore theoretically the possibility 
of such molecular magnetic junctions in which the carbon 
nanocage is a fullerene dimer. Since, among various types 



Author to whom correspondence should be addressed. Electronic 
mail: asaffarz@sfu.ca 



of fullerene dimers, C120 in which the two Cgo molecules 
are connected by a cyclic C4 unit, is the simplest one 
and has the lowest energy as well as interesting physical 
and chemical properties [22] , we have chosen this isomer 
of the Cgo dimer in our calculation. We note that, the 
electronic transport properties of a pure Cgo dimer [23] , 
N and B doped C 60 dimer [24], and Li@C 60 dimer [25] 
have been reported theoretically, however, encapsulation 
of magnetic atoms in fullerene dimers has not been ad- 
dressed in the previous theories. 

We consider an endohedral fullerene Co@Cgo dimer 
contacted via single gold atoms (as in some experimen- 
tally realized single-molecule electronic devices [26, 27]) 
by semi-infinite non-magnetic electrodes, as shown in 
Fig. 1. The electronic density of states, degree of 
spin polarization, spin-polarized charge currents and dif- 
ferential conductances for electron tunneling through 
Au/Co@C6o-dimer/Au junction as a function of bias and 
gate voltage in the dimer ground state, which we find to 
be ferromagnetic, are studied by means of DFT calcula- 
tions of the relaxed geometry and the extended Hiickel 
model of quantum chemistry [26, 28, 29]. The extended 
Hiickel theory (EHT) has been used to study spin trans- 
port in non-magnetic molecules contacted by magnetic 
electrodes [30-32] and in single molecule magnets cou- 
pled to the non-magnetic electrodes [33, 34]. However, 
the EHT does not by itself take account of spin polar- 
ization effects, and for this purpose wc generalized the 
theory to include spin splittings calculated within DFT. 

We carried out ab initio DFT geometry relaxations for 
the Co@Cgo dimer sandwiched between two Au atoms us- 
ing the GAUSSIAN 09 package with the B3PW91 func- 
tional and the Lanl2DZ basis [35, 36]. In the relaxed 
structure (Fig. 1) the Co atoms are at off-center positions 
in the Ceo molecules, with a shortest Co-C bond length 
of 2.09 A. Wc found the molecules to couple to each other 
by parallel bonds with a length of 1.59 A in agreement 
with previous studies [24, 25, 37]. The Au-C and Co- 
Co distances are 2.26 A and 10.49 A, respectively. In 
Fig. 1, the distance between leftmost and rightmost car- 
bon atoms is 16.31 A. Thus our proposed spin-injection 
device has a size less than 2 nanometers. 

The dimer system in the presence of gold electrodes is 
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FIG. 1. (Color online) The relaxed geometry of an endohedral 
fullerene Co@C6o dimer junction. The gray, pink, and yellow 
spheres represent C, Co, and Au atoms, respectively. The 
effect of the gate voltage, Vg, is to shift the molecular orbitals 
relative to the Fermi level of the gold electrodes. 



described by a tight-binding model Hamiltonian derived 
from EHT in a basis of atomic valence orbitals, the 3d, 4s, 
and Ap valence orbitals of the Co atoms and the 2s and 2p 
orbitals of the carbon atoms. These 498 atomic orbitals 
form the basis set for the Co@Ceo dimer Hamiltonian 
that can be written for an electron with spin a as 
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where d^ aa . is the creation operator for an electron in an 
atomic valence orbital Vw with energy eigenvalue e aa . 
"fa/3, a is the spin-dependent matrix element of the ex- 
tended Huckel Hamiltonian between valence orbitals ip aa - 
and ippo-. Since the basis set used in EHT is nonorthog- 
onal, the orbital overlap S a p, a = (ip aa \ip0 a ) that we ob- 
tain from EHT can be non-zero. Therefore, we replace 
7a/3,tr in the Hamiltonian by ja/3,a — eS a p t(7 where e is 
the electron energy [38, 39]. The spin-dependence of 
e aa and 7 Q ^ j0 - results from the spin splittings of the 3d 
and 4s orbitals of the Co atoms that are estimated from 
DFT calculations for major peaks in the spin density of 
states for Co on graphene [40, 41]. We note that the 
spin-dependent electron scattering that arises from these 
exchange related spin-splittings is much stronger than 
that due to spin-orbit coupling. Therefore we neglect 
the latter in the present model. The spin splitting values 
modify the Co 4s and 3d orbitals on-site energies ob- 
tained from the EHT parameters. While the limitations 
of DFT are well known [26] , the present theoretical ap- 
proach has explained [41] the results of scanning tunnel- 
ing spectroscopy experiments [42] on Co atoms adsorbed 
on graphene. Based on the Hamiltonian Eq. (1), the 
total Hamiltonian of the CoQCeo dimer coupled to the 
left (L) and right (R) gold electrodes for an electron with 
spin a can be written as H a = Ho.a + ^L.a + ^R.a- Here, 
S, )lCr (e) = 7V;,erfl l Tj(e)' 7 V7, ( T i s a self-energy that describes the 
coupling of the dimer to the electrode r\ (=L or R), T 7h<7 
is the hopping matrix between electrode r\ and the dimer 
whose matrix elements are given by the EHT parameters. 
For Au atoms in the electrodes we use 6s, 6p x , 6p y , 6p z , 
5d x 2_ y 2 y 5d xy , 5d 32 2_ r 2, 5d zx , 5d zy atomic valence or- 



bitals. Due to the nine valence orbitals for each Au atom, 
we model the electrodes, as in previous work [27, 31- 
34, 43], by nine semi- infinite one-dimensional chains of 
orthogonal atomic orbitals, with one orbital per site and 
periodicity a between orbital sites. The one-dimensional 
chains arc decoupled from each other, i.e., there is no 
hopping between orbitals in different chains. Accord- 
ingly, the surface Green's function matrix for electrode r\ 
is g n = — (1/i) e lka I where k is the electron wave number 
and / is a 9 x 9 unit matrix [44] . We note that, the overlap 
matrix S between the dimer and Au orbitals is included 
in the calculations. In this study, the Co 3d orbitals have 
been shifted in energy so that the average electric charge, 
corresponding to the two spin-up and spin-down states 
of each Co atom in the antiferromagnetic alignment and 
calculated from Eq. 1, is close to the Mulliken charge 
obtained from DFT using the GAUSSIAN 09 package. 
In addition, to locate the position of Fermi energy of 
the gold electrodes between the HOMO (highest occu- 
pied molecular orbitals) and LUMO (lowest unoccupied 
molecular orbital) energies of the CoQCeo dimer, all the 
atomic valence orbitals of the dimer were shifted down in 
energy by 0.4 1 (in units of t = 2.7 cV), relative to their 
energies in the extended Hiickcl model. 

The magnetic moments of the two Co atoms can be 
parallel (ferromagnetic (FM)) or antiparallel (antiferro- 
magnetic (AF)). In the FM case, a spin- up electron trav- 
els from a majority spin state of one Ceo molecule, in- 
duced by its Co atom, to a majority spin state of the 
other Ceo molecule, induced by the other Co atom, while 
spin-down electrons travel from a minority spin to a mi- 
nority spin state. In the AF alignment, electrons be- 
longing to the majority spin state of one molecule, travel 
to the minority spin states of the other molecule and 
visa versa. To find the magnetic ground state configura- 
tion, we estimate the total ground state energy for each 
spin state in the magnetic alignment £ (=FM or AF) 
J2i E^ a where E^ a is the i-th molecular 
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orbital eigenenergy for a spin state and the summation 
is over occupied molecular orbitals. We found the low- 
est total energy to correspond to the spin-up state of 
the FM alignment whose energy is lower than for the 
AF alignment by ~ 5.24 eV. Accordingly, to switch the 
Co@C6o dimer from FM to AF alignment, a strong mag- 
netic field is required. In the absence of a strong field, 
the Co@C6o dimer will remain permanently in its FM 
state. In the present system the strong hybridization of 
the Co 3d orbitals and carbon valence states results in in- 
direct exchange between the Co atoms (mediated by the 
carbon valence orbitals) that stabilizes the FM configura- 
tion that has also been reported in analogous structures 
[45-47]. Therefore, to study the effect of encapsulated Co 
atoms on electronic transport through the dimer bridge, 
only the FM alignment will be considered. 

To calculate the spin states and the degree of spin po- 
larization in this system, we define the spin-dependent 
Green's function of the CoQCeo dimer coupled to the 
gold electrodes as G a (e) = [(e + i5)S a — H^ 1 where 
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FIG. 2. (Color online) Calculated carbon atoms TDOS per 
site vs. electron energy e for the endohedral fullerene Co@C60 
dimer at various gate voltages: (a) Vg = 0.115 t/e, (b) Vg = 
0.0 t/e, (c) Vg = -0.177 t/e. The symbol f (I) corresponds 
to spin-up (spin-down) electrons. The vertical dotted lines 
indicate the position of the Fermi level. All states up to the 
Fermi level are occupied (shaded in the plot). 



5 is a positive infinitesimal. Then the total density of 
states (TDOS) of carbon atoms per site for spin a is 
D a (e) =-j^E a lm((a, a\Ga-S a \a, a)), where the sum 
is over the carbon 2s, 2p x , 2p y , and 2p z orbitals, and 
JV = 120 is the number of carbon atoms in the dimer. 
The degree of density of states spin polarization P is de- 
fined as P = {D t {E F ) - D±(E F ))/(D f (E F ) + D±(E F )), 
where D^(E F )(D±(E F )) represents the TDOS of spin-up 
(spin-down) electrons at the Fermi energy. Note that, the 
spin dependence of electronic states at the sites of carbon 
atoms originates from hybridization between carbon or- 
bitals and Co 4s and 3d orbitals. To see the effect of this 
hybridization on the electronic states, we show in Fig. 2 
the TDOS as a function of electron energy for represen- 
tative gate voltages Vg = +0.115, 0.0, and —0.177 t/e. 
In our calculations, the extended Hiickcl model param- 
eters describing both the on-site and intersite Hamilto- 
nian matrix elements are recalculated for each value of 
the gate voltage. The dependence of the intersite Hamil- 
tonian matrix elements on the gate voltage (see Ref. 30, 
Appendix A) implies that hybridization between the or- 
bitals on different atoms depends on the gate voltage in 
the present theory. However, since the same gate poten- 
tial is applied to the Ceo molecules and Co atoms, the 
molecular orbital features shift nearly rigidly under gat- 
ing. The chosen positive and negative gate voltage values 
in Fig. 2 correspond to the first peaks (in the absolute 
value of P in the positive and negative region of gate 
voltage) crossing the Fermi level (see the arrows in Fig. 
3). At zero gate voltage the Fermi level is located in the 
HOMO-LUMO gap where no molecular orbital is avail- 
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FIG. 3. (Color online) Degree of density of sates spin polar- 
ization P at the Fermi energy as a function of gate voltage. 
The arrows show the positions of the first maxima of spin 
polarization with varying gate voltage. 



able for electrons to tunnel through the dimer (see Fig. 
2(b)), and hence little or no spin polarization is expected. 
By applying a gate voltage, the encapsulated Co atoms 
are ionized due to the gate-induced shift in energy of the 
molecular orbitals relative to the gold electrode's Fermi 
energy. The effects of the positive and the negative gate 
voltage, have been shown in Fig. 2(a) and 2(c), respec- 
tively. Comparing Fig. 2(b) with Fig. 2(a) and 2(c), it 
is clear that a positive (negative) gate voltage shifts the 
molecular orbitals to a higher (lower) energy and for the 
value of Vg = +0.115 t/e at which the spin down states 
are conductive, there is a significant difference between 
the two spin populations which can cause a strong spin 
polarization at the Fermi level as shown in Fig. 3. On the 
other hand, in the case of Vg = —0.177 t/e the molecular 
orbitals shift to a lower energy and, as shown in Fig. 3, a 
peak in the spin polarization appears. The results shown 
in Fig. 3 for P versus Vg indicate that the spin injec- 
tion can be controlled by changing the gate voltage and 
positive or negative values for the spin polarization may 
be obtained, depending on the value and the sign of gate 
potential. The magnitude of P can exceed 90% for cer- 
tain values of gate voltage which is considerably higher 
than that calculated for Co-nanowirc@CNT (66%) [5] 
and those measured for Co-cluster/Cgo mixtures (70%) 
[7], Co-nanoparticle/Co-C6o compounds (40%) [11] and 
Co crystals (42%) [48]. 

To investigate coherent charge transport through the 
Co@C6o dimer, we make use of the Landaucr-Buttikcr 
formula based on the nonequilibrium Green's function 
method [44] for which the spin-polarized charge current 
at a constant bias voltage, Vb, is calculated as 

/oo 
T„(e,V b )[f(e-n L )-f(e-n R )]de, (2) 

where /(e) is the Fermi function, hl,b, = £f ± feVf, 
are the chemical potentials of the electrodes, T a = 
Tr[r.L ]0 .G .r-R ]0 .Gj.] is the spin-dependent transmission 
function for electron tunneling through the dimer, and 
r r)]Cr = — 2Im(E r?i(T ) is the coupling matrix between the 
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FIG. 4. (Color online) Calculated spin-dependent electric 
currents [(a), (b)], current spin polarization Pi (c), and 
differential conductances [(d), (e)] vs bias voltage for the 
Au/Co@C6o-dimer/Au junction, shown in Fig. 1, at various 
Vg values. 



dimer and the electrode rj. 

The calculated current-voltage characteristics, I-V, 
and the differential conductance spectra, dl/dV, of the 
Au/Co@C6o-dimer/Au junction as functions of bias volt- 
age Vb at Vg = 0.0, 0.1, 0.2, and 0.3 t/e and at temper- 
ature T = 4.2 K for spin-up and spin-down electrons are 
shown in Fig. 4. Due to the structural symmetry in the 
system, we obtained a symmetrical behavior in the I-V 
curves with respect to Vb = 0.0 t/e. Therefore, the results 
for negative bias voltages are not shown in Fig. 4(a) and 
4(b). The spin up and spin down charge currents show a 
step-like behavior due to the entry of molecular orbitals 
into the bias window. The sharpness of the steps depends 
on the hybridization of the dimer electronic states shown 
in Fig. 2 to the gold electrode orbitals. The broadening 
of the quantized states in the scattering region causes 
smooth steps in the I-V curves for both spin directions. 
We see that the I-V characteristics of the Co@Cgo dimer 
reveal a conducting or semiconducting behavior depend- 
ing on the value of gate and bias voltages. This electronic 
property is strongly spin dependent due to the spin split- 
ting of Co 3d orbitals [40, 41] around the Fermi energy of 
gold electrodes. This is consistent with the work of Xic 
et at [17] who found that hybridization between Co 3d 
and C 2p orbitals for Co nanowires encapsulated in CNTs 
is strong and results in spin polarization that increases 
with increasing thickness of the nanowire. 

At V G = 0.0, 0.1, and 0.2 t/e, a threshold voltage that 
is different for spin-up and spin-down electrons, is needed 
to generate a charge current through the device. That is, 
at low gate voltages for low bias voltages the device is in 



an off state, while at Vg = 0.3 t/e, the device is turned 
on and the charge current increases linearly, indicating 
Olimic behavior at low biases. Such a behavior is rea- 
sonable because in this situation, the gate voltage shifts 
the molecular orbitals within the bias window and more 
channels for tunneling electrons through the Co@C6o 
dimer bridge become accessible. Furthermore, by com- 
paring the charge currents in Fig. 4(a) and 4(b) at each 
bias voltage, one can readily confirm that the tunneling 
channels for spin-up and spin-down electrons are well sep- 
arated due to the presence of Co atoms inside the dimer. 
Thus, our model predicts gate and bias voltage-controlled 
switching between off and on states and also between 
spin-up and spin-down electric currents. In addition, the 
results shown in Fig. 4(a) and 4(b) reveal that the encap- 
sulation of Co atoms in fullcrcnc cages is an efficient way 
for injecting spin-polarized electrons into a nonmagnetic 
conductor by electric current. To estimate the spin injec- 
tion efficiency in our model, we calculate the current spin 
polarization defined as Pj = (If — ij,)/(/f + Ij.) for each 
bias and gate voltage. The results for Pj as a function of 
V, are shown in the Fig. 4(c). We clearly see that very 
high values for current spin polarization (~ 99%) can be 
obtained in this hybrid structure. These values depend 
strongly on the bias and gate voltages and predict that 
highly efficient voltage-controlled spin injection can be 
achieved using the Co@C6o dimer. 

The differential conductances shown in Fig. 4(d) and 
(e) clarify the role of gate and bias voltages in control- 
ling the spin-polarized transport through the system. To 
show the conductance gap, i.e., the voltage width of 
the zero-conductance region, we have plotted the dl/dV 
spectra for both positive and negative bias voltages. In 
Fig. 4(d) and (e) the size of the conductance gap de- 
creases and the system gradually moves into a conducting 
state with increasing Vq . The symmetry of the differen- 
tial conductances about Vb — confirms the existence of 
symmetry in the I-V curves of Fig. 4(a) and (b). The 
sharp conductance peaks correspond to the step-like fea- 
tures seen in the I-V curves and signal the opening of 
new conducting channels through the dimer. 

In conclusion, we have shown, based on ah initio and 
semi-empirical calculations, that very high degrees of 
spin polarization can be achieved for electrons tunneling 
through cndohcdral fullerene Co@Cgo dimers in their fer- 
romagnetic ground state. Our results for spin-dependent 
I-V curves and dl/dV spectra, calculated by means of 
Landaucr-Buttiker theory, show that these hybrid struc- 
tures can be utilized for highly efficient spin injection 
into a nonmagnetic conductor by adjusting gate and/or 
bias voltages. This electrostatic activation should lend 
itself to integration of such nano spin injectors with con- 
ventional nano electronic devices. Also, since magnetic 
electrodes are not required to achieve spin injection in 
this system, the size of the device required for spin in- 
jection is that of the Co@Cgo dimer itself, i.e., less than 
2 nanometers. Therefore as well as being of fundamen- 
tal interest, systematic experimental studies of encapsu- 



latcd magnetic atoms in the fullercnc dimcr based on 
these findings may be relevant for nanoelectronic and 



spintronic applications. 
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